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INTRODUCTION 
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I. 

Optimization  problems  have  been  investigated  the  past  several  years 
by  many  researchers.  Uzkov[l]  first  found  the  solution  to  the  directive 
gain  maximization  by  employing  linear  transformations.  Rloch,  et  al.[2], 
Morgan[3]  and  Uzsoky  and  Solymar[4]  gave  some  applications  for  arrays 
with  arbitrary  spacings.  Tai[5]  considered  the  problem  of  achieving 
maximum  gain  in  uniform  linear  arrays  of  dipoles  and  gave  some  graphs  on 
the  subject.  The  optimization  of  SNR  with  constraints  on  the  Q-factor 
has  been  obtained  by  Lo,  et  a 1 . [6] . 

Cheng[7]  studied  the  maximum  directive  gain  by  the  eigenvalue 
method.  Sunzgiri  and  Butler[8]  have  used  the  eigenvalue  method  to  find 
the  solution  for  the  maximum  directive  gain  with  constraints  on  the  re- 
sulting sidelobes.  Matrix  methods  were  applied  by  Strait  and  Kuo[9]  for 
constrained  optimization  of  various  performance  indices  of  arrays  of 
straight,  parallel  thin-wire  antennas. 

Sa ha  1 os [ 10]  gave  a solution  for  maximization  of  the  directive  gain 
of  arrays  with  arbitrarily  oriented  short  dipoles.  Sarkar  and  Strait[14] 
gave  some  useful  formulas  for  optimizing  certain  performance  indices  of 
arbitrarily  oriented  arrays  of  wire  antennas.  All  these  formulas  were 
found  assuming  known  polarization  in  the  specific  direction  for  which  the 
maximum  of  index  is  desired,  and  their  study  is  for  optimization  problems 
with  only  one  kind  of  constraint  (i.e.,  nulls  or  sidelobes  or  other 
indices  but  not  more  than  one  simultaneously). 

This  paper  presents  some  unified  formulations  for  optimizing  any 
performance  index  subject  to  constraints  on  other  indices,  sidelobe 
levels  and  pattern  nulls  of  arrays  with  nonparallel  thin  wire  antennas. 
The  formulations  presented  here  are  not  restricted  to  optimizations  with 
only  one  constraint  at  a time  nor  do  they  assume  the  polarization  is 
known . 


II.  MOMENT  METHOD  FORMULATION 

The  moment  method  in  electromagnetics  is  very  useful  because  it 
transforms  problems  which  are  formulated  in  terms  of  integral  equations 
into  systems  of  simultaneous  linear  equations.  Let  a wire  structure  be 
composed  of  a number  of  straight  segments.  If  we  define  a right  handed 
orthogonal  coordinate  system  (n,t,f)  at  each  point  of  the  wire's  cylindri- 
cal surface,  we  can  find  [11]  that: 

- L 1(0  (E"'-ZsHn,)d*  = Vm 
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Equation  (1)  is  the  reaction  integral  equation  developed  by  Richmond  and 
is  true  for  electrically  thin  wires. 

!(■)  is  the  total  current,  t is  a metric  coordinate  measuring  position 
along  the  wire  axis,  is  the  surface  impedance  for  exterior  excitation 
and  L is  the  wire  length.  The  only  unknown  in  Eq.  (1)  is  the  current 
!(>)•  The  other  quantities  are  given  as  follows: 
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are  the  impressed  currents  by  which  the  wire  structure  is  excited. 

(Em,Hrn)  is  the  field  of  an  electric  test  source  located  in  the  interior 
of  the  wire  surface  and  radiating  in  free  space. 

The  next  step  is  to  expand  the  current  on  the  wire  in  a finite 
series  of  the  form: 


1(0  = l InFn(0 


(5) 


By  Eqs.  (1),  (5)  we  can  write  a system  of  simultaneous  linear  algebraic 
equations 


[Z]  (I)  = (V) 


(6) 


where  the  elements  Zmn  of  the  matrix  [Z]  are  of  the  form 


Zmn  ■ • j F„(0(E^-Zs^)d. 


(7) 
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Solving  Eq.  (6)  for  the  current  we  have 

(I)  = [Z  r](V)  = [Y]  (V ) . (8) 


In  an  antenna  array  with  N elements,  we  assume  that  each  element 
is  fed  at  only  one  port.  rrom  Nj  elements  of  the  column  vector  (V), 

(N|  N-Np,  where  Np  is  the  number  of  the  expansion  modes),  only  N are 
non-zero.  Thus,  we  can  find  a matrix  equation  which  relates  the  currents 
on  the  main  ports  with  the  corresponding  voltages  of  the  form: 


(!')  = [Yr](V') 


(9) 


The  matrix  [Yr]  contains  the  first  N rows  and  columns  of  the  matrix  [Y] 
which  are  associated  with  the  non-zero  entries  in  (V). 


III.  FORMULATION  OF  THE  PERFORMANCE  INDICES 
A.  Antenna  Far  Fi eld 

In  finding  the  solution  of  an  EM  optimization  problem,  it  is  first 
necessary  to  define  the  total  field  that  is  involved. 

In  the  present  situation  we  are  interested  in  the  total  electric 
field  radiated  by  an  antenna. 


E ( r , <(> , e ) = E0(r,'t>»o)e  + E^r.^.e)* 


(10) 


The  components  E0(r,<j>,o)  and  E^( r,$,fl) 


can  be  expressed  as 


Eg(r  ,<l,0) 


N1 


InEnoU.e) 


N 


7 l 

r i=l 


N1 

l,  YniViEne(^0) 
n=l 


(ID 


1 N1 

E.(r^.e)  = 7 l InEnt)(^9) 
♦ n=l 


1 N 

L y 

r 


i = l 


N1 

l YniviW*’°) 
n = l 


(12) 


where  Yn^  are  the  elements  of  the  matrix 
We  define  two  row  vectors  [B|]  and  [B^], 
each  is  given  by 


[Y],  which  is  given  in  Eq.  (8). 
such  that  the  1th  component  of 
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(13) 


Bli  = l Vni[n  (4>»0 ) 

n ni  ne 


!2i  * I.  Vr,iEn,<»'s> 
n=l 


By  the  above  Eqs.  (10)-(14)  we  can  write  the  far  zone  field  as: 


> i N N 

E(r,^,-)  = I l B -j  • V • e + l B2.V  * 

r i = l i=l  n 1 


^ ( r » f » 0 ) = l C [ B -j  ][V 1 ]o  + [B2][V]i] 


B.  Power  Gain 

The  power  gain  of  an  antenna  is  defined  by 

G = ^(radiation  intensity  for  the  specified  direction) 
Power  input  to  the  array 

In  relation  to  the  electric  field  (17)  can  be  written 


GU.e)  - rZ/30 

pi  n 


The  power  input  can  be  expressed  as: 


Pin  = Re 


(j,  ■ [»•].  ^ tv] 


where  -v  signifies  the  transpose  and  the  * signifies  the  complex  conjugate. 
Equation  (18)  becomes 


- ’C'r="r'  - -wr' 


(20) 


CV‘]*C[b'i]*[b13+[S2]*[b2J][v*3 

G ( » 0 ) - 7" 

1 5[V 1 ]*[[Yr]+[Yr]*][V] 

C . Directive  Gain  and  Direct i v vty 
The  directive  gain  is  defined  by: 

n _ 4n(radiation  intensity  for  the  specified  di recti onj 
radiation  power 

Equation  (21)  may  be  expressed  as 

d ( j , o ) = (22) 

Pr 

and  Pr  = P^n  - Pq,  where  Pq  is  the  total  power  dissipated. 

It  can  be  shown  that [12] 

Pr  = Re[[V']*[[Yr]-[Y]][V']]  . (23) 


Complete  details  for  the  elements  of  [Y]  are  available  in  Reference  [12]. 
The  maximum  value  of  the  D(d>,f>)  is  the  directivity. 

D.  Other  Indices  and  factors 

Efficiency  indices  and  the  sensitivity  factors  are  very  important 
parameters  describing  the  performance  of  an  array.  Efficiency  indices 
can  be  defined  in  several  ways.  Two  of  these  are: 


_ radiation  intensity  in  the  direction  of  max,  radiation  /£4) 
sum  of  the  excitation  voltage  magnitudes  squared 


s = raAi At i 0JL  intensity  j n the  direction  of  max,  radiation 
^1  sum  of  the  feed  port  current  magnitudes  squared 


(25) 


Equations  (24)  and  (25)  can  be  written: 

s [v  ] *[b1]-«-[b2]*cb23]  [V] 

1 20^[V ’ ]*[V  * ] 


(26) 
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(27) 


[^' ]*[[B1 ]*[B1 ]>[B2]*[B2][V ] 
1 20n ' ]*[yr]*[Yr][V'] 


Correspondingly,  sensitivity  factors  can  be  defined  as: 

i 1 . . 1 

k = s andkrsj  • 

Another  index  is  the  Q-factor  that  relates  gain  to  efficiency  index  or 
sensitivity  factor.  These  are  defined  as: 


Q = G/S  = G-K 

(28) 

and 

Q]  = G/S-)  = G K-] 

(29) 

IV. 

ARRAY  OPTIMIZATION 

A.  Gain-Directive  Gain  and  Efficiency 
Index  Maximization 

As  we  have  seen  the  power  gain  is 

GU,e)  = - ! - - (30) 

[V  ■ ]*[  M2][V] 


where 


[M2]  = 15  [ [Yr]+[Yr]*]  . (31) 

Taking  the  first  variable  of  0,63,  equal  to  zero  we  can  find  the  [V‘] 
which  maximize  G.  Thus, 
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AG  = [ ' V ' ] * 


[B,]< 


[B,][V] 


[V,]*[M2][V] 

[^’ 3*Crft13*[B13f[g23*CB233Cv] 

[V]*[M  2][V']I? 

. ★ 

•!  I [«V]  = 0 


[B?]’ 


CB 2][V‘] 
[V']*[M2][V] 


[M2][V] 


(32) 


The  quantity  { * is  the  conjugate  transpose  vector  of  the  first  { f. 

By  Eq.  (32)  we  have 

f 


Re  -j  MO 


(33) 


and  so  because  we  are  interested  only  in  relative  values  of  voltages, 
we  have: 


[v]  = [M2r]  [[&,]*  + *[b2]*] 


(34) 


where 


[B2][V] 

a = = a value  dePendent  on  the  polarization  in  the 

*-  1 J L v J direction  of  radiation. 


With  the  same  procedure  we  can  find  the  feed  voltages  for  maximization 
of  the  directive  gain  or  any  efficiency  index.  The  difference  will  be 
in  the  matrix  [Mp],  which  for  the  above  indices  we  can  find  in  Eqs.  (23), 
(26),  and  (27). 


B . Optimization  Subject  to  Constraints 
on  Radiation  Pattern  Nulls 

Any  index  such  as  for  gain,  directive  gain,  or  efficiency  can  be 
written  in  the  general  form: 
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(35) 


Ctf,3*CC81]*[B13^EB2]*[B23][v»] 

[V']*[M2][V] 


Suppose  that  it  is  desired  to  determine  feed  voltages  that  will 
provide  pattern  nulls  in  (P-1)  directions  and  maximum  , subject  to  this 
constraint  in  the  pth  given  direction.  If  Eq.  (16)  is  written  for  each 
of  the  P directions  we  will  have: 


'to 

(fl  “>+f2  4) 

K 11  ir; 

0 

^f1  • ff?  * ) 
'12  T12' ' 

0 

= 

* 

0 

V 


2 


Equation  (36)  may  be  modified  to: 


(36) 


0 

0 


'"f1  f1 
11  12 

f2  f2 
11  12 


IN 

1h 


(37) 


or 


[0] 


where 


[V|]  column  vector  with  N-2(P-1)  elements 
[ ] column  vector  with  2(P-1)  elements 


(38) 


I 


\ 


- 
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[B^j]  matrix  with  [N-2 (P-1 ) ] x 2(P-1)  elements 
[Bpp]  matrix  with  2 ( P- 1 ) x 2(p-l)  elements 


[E0J  = [B11  j B12^  v, 


From  Eq.  (38)  we  have: 


[v$]  = - [B22J-  [B21][V|]  = [D][V| ] 


[E0]  = [[Bn]+[B12][D]]  [V]] 


We  easily  can  see  that  EQ  is  of  the  same  form  as  Eq.  (16). 


t0  - [[r1]H+[r2j;][V'] 


where: 


[r^  = [[8nH^2][D]]-o 
[ 1 2 = ft  Bl  1 ] + [^i  2][D]  ] ' 


From  Eqs.  (40)  and  (42)  the  index  y becomes 


C^i]*[[r1D*[r1]+Cr2D*[r2]][v^] 


[Vi  ]*[Q][V-j  ] 


where 


[Q]  = [U  i -b*][M2] 


([U]  is  a unitary  matrix). 
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Now  the  maximization  of  index  r under  constraints  of  radiation 
nulls  is  achieved  by  maximization  of  i.e.,  Eq.  (44). 

With  the  same  procedure  as  before  we  can  find: 


[V]]  = [Q]_1[[:'i3**a[:2]*] 


and 


'max 


[r1][Q]'1[[r1]*+^r2]*] 


(46) 


(47) 


C . Optimization  Subject  to  Constraints 
on  Side-lobe  Level 


Next,  consider  the  problem  of  finding  [V1]  that  will  maximize  the 
index  y of  Eq.  (35)  in  a given  direction  subject  to  a constraint  on  the 
sidelobe  level.  Suppose  that  we  have  p known  sidelobes  in  p directions. 
That  means: 


-► 

El 

-X. 

Al 

E 2 

k2 

a2 

= 

E 9+ 
oe 

• 

_ 

_kP. 

_AP_ 

E-off  [K]Eq(iG+[-,']Eoi£4| 


(48) 


where  Eo0  and  E0.  are  the  components  of  the  electric  field  £0  in  the 
direction  that  maximizes  the  index  y. 


From  Eq.  (48)  we  can  take 


[K]Eo0 

[ A ] E 


0*  j 


[B2i  i b22] 


V] 

V* 


(49) 


and 


[v$]  = [b22]_1 


WQI) 

Wo, 


- [B21][V]] 


(53) 
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From  Eq.  (39)  we  have 


[£o]  = [^11  1 ^1?] 


So  Eqs.  (50)  and  (51)  give: 


[f0l  = [Bn]-[B12J[B22]-1[B2)]  [V- ]+[B12][B22]-1  0 | (52) 

^ ^ ^ r\ 


r [DfJCV  3^[tt 1 2 v22^ 


Eo*  = ] + [”i2  71 22 ^ 


F[K]Er 


[D../^]  = [ [B]  ] ]-[B] 1 ] ' "/* 


[ 71 1 2 1 71  -1  = [^12^^22^  ' 


The  system  Eq.  (53)  gives 


E0  = [[e1]o+[;:2]^][V1] 


£ ii 

°22l  J 2 ] L • J 

a22aira21al? 

i - 

-a12[D0]+a1,[D} 

a2?aira21a12 

'll 


‘12 


i ■ [ n ] 2 ] r k ] 

- c-;2][a] 

a2]  ' " 
a 22  1 1 ' 

With  the  help  of  Eqs.  (50)  and  (53),  [V^]  takes  the  fonn: 

r : 


•1 


[K][  t] 

I 1 

LCa][j:2] 


- [B21] 


[Vj]  = [D][Vj] 


(56) 


(57) 


The  performance  index  with  the  help  of  Eqs.  (55)  and  (57)  becomes 


[Vj]*[Q][vj] 

and  from  the  procedure  as  before,  the  maximum  y is  obtained  when 
[V]]  = [Q]_1[[ : 1]*+a[>:2]*] 

where  [Q]  is  given  by  Eq.  (45). 

The  above  procedure  gives  only  the  maximum  index  y when  we  know  p 
values  of  the  electric  field.  For  the  sidelobe  levels  to  be  the  desired 
values  we  must  use  an  iterative  procedure.  So,  we  start  from  a known 
radiation  pattern  which  has  the  desired  sidelobes  and  use  the  values  of 
this  pattern  in  the  above  procedure.  With  this  we  can  find  feed  voltages 
[Vj]  and  the  level  and  directions  of  sidelobes.  These  sidelobes  are  com- 
pared with  desired  results.  The  new  directions  of  sidelobes  are  used  to 


(58) 


(59) 
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tin.)  new  teed  voltages  [Vj],  etc.  The  iterative  procedure  is  continued 
until  all  sidelobes  will  be  of  the  desired  values. 


D.  Example  of  a Radiation  Pattern 
Useful  as  the  Starting  Point 
for  the  Iterative  Procedure 

Suppose  that  we  want  to  find  an  antenna  which  has 

i . j ict-vg.^/z)!  . k 

IE*(i/2.»/2)! 

(ii)  ~~  lobe  ratio  > ?- 

' Side  4 


We  design  a linear  uniform  array  with  infinitesimally  small 


E(*,n/2)  -2  ^ I ■jCos,t ' cos^~ -J-  dcos^J  • e + 

N . v _ 

+ 2 sint  l I^sino1  cos( dcos^j  ^ 
i=l  V * ' 

This  electric  field  is  that  of  a Dol ph-Chebyshev  array 

we  put 
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dipoles 


(60) 
So  if 


(61) 


N 

Ii  «-0S  ’ = T^n_  , (ax ) 

N I 

^ 1 i S1n  T?i-l(x)  = T2N_i (bx) 

x "(cos^l  cost)) 

> 

we  obtain 

\ = T2N-1 (d)  . . 

Rb  T2N-1^ 

9 

For  the  maximum  value  of  sidelobe  to  be  2, 

R2  + R2  = R2 
a b 

From  Eqs . (62)  and  (63)  we  have: 


R = R 

Ka  -7=5= 

Vk2+T 


From  Ra  and  R^  we  can  define  a and  b as: 


a 


J/2N-1 


-M-’J 


[1/2N-1 


We  can  find,  from  the  Dolph  procedure,  the  Ijcoso1  and  I, 
naturally  the  I-j.n1. 


More  analytical  ly,  if  we  want  to  have  an  array  with 
equal  distance  a/2  and  k = /3,  R=40,  from  Eq.  (60) 

t 


(62) 

(63) 

(64) 

. • (65) 

sinn1  and 
8-dipoles  of 
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Ra  = 36.64  , Rb 
and  by  Eq.  (65) 

a = 1.2  , b = 1.15 


The  results  are  shown  in  Table  I. 

TABLE  I 


I|Sino^  = 3.1 

I-|Cos'^  = 9.85 

l2Sino2  = 2.6 

Ipcose^  = 3.84 

I3sino3  =1.7 

I3coso^  = 0.595 

I^sino4  =1.0 

I^coso^  = 0.28 

TABLE  II 


I = 30.6 

n1  = 5°. 8 

I2  = 4.64 

o2  = 34°. 1 

I3  = 1.80 

03  = 70°. 7 

o 

ii 



CD 

II 

O 

An  appropriate  design  now  would  be  to  use  an  array  of  wires  with  the  side- 
lobes  in  the  same  directions  as  those  of  the  array  of  infinitesimal  dipoles. 


E . Optimization  of  One  Performance 
Index  Subject  to  a Constraint 
on  AnojJier 

Suppose  that  we  want  to  maximize  the  index  y1  subject  to  a constraint 
on  another  index  Y2*  The  two  indices  are  of  the  general  form 
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(66) 


[V']*[  [Bi3*[  B1h[B2]*[B233  [ Vj 
Y1  = ['<f,]*[M2][V'] 


[V 1 ]*[M3] [V 1 ] 
'2  = [V']*[M4XV;] 


(67) 


By  the  Lagrange  method  for  maximizing  y-j  when  y2  6as  a given  value,  we 
must  find  voltage  [V1]  which  maximizes 

[V,]*[[Bi]*[B1]^[B2]*[B2]][Vi] 

[V,]*[M2][V] 


+ 


[v*  3*[m3][v 1 ] 
Cv' ]*[m4][v  ] 


(68) 


With  the  same  procedure,  as  in  IV-A,  we  can  find  that: 

[V]  = [K]'1  [[B1]*+g[B2]*]  (69) 

where 


[K]  = [[M2]+P[Y2[M2]-[M3]]]  . (70) 

The  unknown  value  p we  can  find  with  the  help  of  the  constraint 
Eq.  (67).  By  a method  similar  to  that  of  Winkler  and  Schwartz[13]  the 
substitution  of  [V1]  into  Eg.  (67)  can  give  P from  the  real  eigenvalues 
a of 


[G][X]  = a[X]  (71) 

where 

1 

a = p 
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One  of  the  real  eigenvalues  gives  the  optimum  y,  subject  to  the 
correct  constraint 


F . Op tim i z a tion  of  One  P er f orma nee 
Inde_x_  Subject  to  Constraint 
on  Other  Indices 

Next,  we  consider  the  more  general  problem  of  maximization  of  the 
index  yj  subject  to  constraints  on  other  indices  Y^ » r 3 , • • • , y^ • Index  y, 
is  of  the  same  form  as  before,  and  the  other  indices  are  given  as  follows 

Yi a fn*fMi+1][v]  • {72) 

By  the  method  of  Lagrange  we  can  find  again  that  [Appendix  I] 

[V1]  = [K]-1[[B1]*+u[B2]*]  (73) 


where: 


[K]  = [M2]+P1[y2[M4]-[M3]]  + 

+ PnCynC^^j-CMnDD  • (74) 


The  values  Pi  can  be  found  by  substituting  Eq.  (73)  into  Eq.  (72),  the 
constraint  equation. 


G.  Optimization  of  One  Performance 
I ndex  Subject  to  Constraint  on' 

Other  Indices,  Pattern  Nulls 
ancfTT  cfelob  e Lev  el  s 

As  we  have  seen  in  previous  sections,  pattern  nulls  or  sidelobe 
levels  give  a linear  relation  between  the  voltages 


[V$]  = [D][Vj ] 


(75) 


So,  the  index 


rl 


becomes  as  follows: 


(76) 


[vjMEr^n^MrpnrglHvj] 

M Ev[}*[Q23[vi] 

where 

EQ2]  = 

Equation  (75)  transforms  all  indices  Y2’Y3’"YN  as 

E^MQiJEvj] 


. -V1 

1 

cz 
1 

u ; d* 

[M2] 

— 

• 

D 

Yi  = [vj]*[Qi+1][Vj] 


(77) 


(78) 


where 


[Qi]  - 


1 'V/ 

U ! D* 


[Mi] 


(79) 


Our  problem  of  optimization  of  one  index  subject  to  constraints  on 
other  indices,  pattern  nulls  and  sidelobe  levels  is  achieved  with  the 
help  of  Eqs.  (75),  (76)  and  (78).  We  can  find  the  vector  [ V 4 ] by  the 
help  of  Eqs.  (73)  and  (74)  which  will  be  written  as: 


[Vp  = [K]-1  [[f1]*+u[r2]*]  (80) 


and 


[K]  - [Q2]-P,[Y2tQ4]-[Q3]]  * • ■ • ■ 

WWii-¥  • (fi,> 

Equations  (80)  and  (81)  give  a solution  to  the  more  general  optimization 
problems. 
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V. 


POLARIZATION 


One  of  the  main  problems  in  our  study  is  the  polarization  in  the 
specific  direction  for  which  the  maximum  of  one  index  is  desired.  Sarkar 
and  S t ra i t [ 1 4 ] assume  that  this  is  known.  The  truth  is  that  only  arrays 
with  parallel  wire  antennas  give  us  a known  value  of  polarization.  An 
array  with  nonparallel  wires  has  a specific  value  of  polarization  which 
we  must  find.  Equation  (80)  gives  a general  formula  for  the  voltages 
[Vj].  The  number  A dependent  on  the  polarization  is  given  by: 


r2][Vj] 

[TTTFjl 


From  Eqs.  (80)  and  (82)  we  can  find  that  A is  the  solution  of  the 
equation: 


a2[t1][k] 


+ >{[  - [i 2][K]_1Cr2]*>  + 

+ [r2][K]  [?,]*  = o 


(83) 


When  we  have  only  linear  relations  between  the  voltages,  i.e.,  constraints 
on  nulls  and  sidelobes,  Eq.  (83)  gives  directly  the  value  of  . Con- 
straints of  other  indices  give  nonlinear  relations  between  the  voltages 
and  in  this  case  an  iterative  method  can  give  us  the  a.  At  first  we 
suppose  that  A is  known  and  given  by  Eq.  (83)  where  the  matrix  [K]  is 
equal  to  [Qo]  from  Eq.  (77).  Using  this  value  of  A we  can  find  matrix 
[ K]  by  whicn  Eq.  (83)  gives  a new  The  new  A is  used  to  find  a new  [K], 
etc.  The  iterative  procedure  is  continued  until  the  new  \ is  the  same  as 
the  previously  found  one.  Numerical  results  obtained  by  computer  showed 
that  there  are  two  equivalent  values  of  \ which  give  the  same  maximum 
index.  Each  one  of  these  is  the  negative  inverse  of  the  other. 


VI.  EXAMPLES 

1.  Consider  the  design  of  6 and  10-element  linear  arrays  of  parallel, 
half-wavelength,  centerfed  wires  to  provide  optimum  indices  in  the  broad- 
side direction  ( :o=h/2,u0=tt/?) . The  wires  are  assumed  all  of  radius  0.001\. 
Figure  1 shows  the  normalized  E°(*,i72)  and  E*( 4. ,ti/2 ) fields  of  a linear 
array  with  6 elements  in  equal  distance  A/2  and  direction  <f>1=0,i'»1=n/4. 

The  maximum  gain  is  G= 10.97  subject  to  the  constraint  that  nulls  are 
required  in  both  the  direction  ($=66  ,o=90°)  and  the  direction  ( g>=48  , 
0=90°).  A 10-element  linear  array  with  the  same  geometry  as  above  has 
maximum  gain  16.00  subject  to  the  constraints  that:  (i)  nulls  are  required 
in  the  directions  ($=58. 6" ,0=90'  ) , ($=34'’, 0=90°) , (ii)  the  E"(”/2,i)  side- 
lobe  levels  are  1/5  the  level  of  the  main  beam  and  (iii)  the  efficiency 
index  S is  48.0.  Figure  2 shows  the  normalized  E°U,n/2)  and  E^(i(>,n/2) 
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fields.  A 6-element  linear  array  with  the  elements  in  directions  f'  0, 

, = (i-lW10  and  equal  distance  '/2,  has  optimum  broadside  polarization 
P=1.720  or  P=-.579  subject  to  the  constraint  that  a null  is  required  in 
the  direction  t=50o,e=90°.  The  maximum  gain  is  10.24  and  in  Fig.  3 we 
can  see  the  normalized  E‘'(.t,n/2)  and  E^(<t>,n/2)  fields. 

The  same  6 element  array  is  considered  subject  to  constraints  that: 

(i)  a null  is  in  the  direction  4=55° ,0=90° , (ii)  the  E',('t>,n/2)  sidelobe 
levels  are  less  than  1/10  and,  (iii)  the  efficiency  index  is  S=65  with 
optimum  polarization  P= 1 .47  or  -0.68,  and  maximum  gain  8.89.  Figure  4 
shows  the  E "(t,n/2)  and  ET(;,ti/2)  normalized  fields. 

If  nulls  were  required  in  the  directions  t=60n ,<»  = 90'  for  the  E and 
$=40°, 0=90°  for  the  E*  field,  then  the  maximum  gain  goes  to  11.2.  In 
Fig.  5 we  can  see  the  fields  in  o=tt/2  plane. for  a 10-element  linear  array 
with  the  elenents  in  directions  t1=0  and  01  = ( i -1  )n/18  and  equal  distance 
>/2  subject  to  constraints  that:  (i)  there  are  nulls  in  directions  :=75°, 

u=90°  for  the  EJ  and  ^=60° ,e=90°  for  the  ET  field,  (ii)  the  sidelobe 
levels  of  E"($,ti/2)  field  are  less  than  (1/10)  and,  (iii)  the  efficiency 
index  is  S=71.0  with  optimum  polarization  P= 1 . 542 , or  -.648,  and  maximum 
gain  G=20.64. 

The  same  array  is  considered  with  (i)  nulls  in  directions  *=70° ,0=90° 
for  EfJ  and  *=35n ,0=90°  for  E*  field,  (ii)  gain  G=20  and,  (iii)  E“(:,n/2), 
sidelobe  levels  less  than  1/12.5,  has  optimum  polarization  P=1.796,  or 
-.557,  and  maximum  efficiency  index  S=94.38.  Figures  6 and  7 show  the 
E0U,"/2)  and  E4, ( t /2 > normalized  fields  for  the  above  two  cases. 


2.  Directivity  versus  the  interelement  distance  d/A  for  a 10-elem<-nt 
linear  array  is  shown  in  Figs.  8a  and  b.  Elements  have  the  directions 
1 =0 »o‘>  = (i-l  )^/18.  In  Fig.  8a  we  can  see  the  maximum  directivity  under 
constraints  on  nulls  in  directions  ( <t>=75°, 0=90°) , (t=40°, 0=90°)  for  En 
and  ( 4>=60°, 0=90° ) for  the  E*  field.  Figure  8b  shows  the  maximum  direc- 
tivity under  constraints  on  sidelobe  levels  as  follows: 

(i)  114®  lobe  = for  the  E9  field. 


(ii)  ^4^  lobe  = is-F  for  the  E4  field. 

' 'main  12.5 

The  corresponding  optimum  polarization  for  both  cases  is  represented  in 
Figs.  9a  and  b.  The  maximum  directivity  subject  to  constraints  on  nulls 
for  two  different  10-element  arrays  we  have  on  Fig.  10.  Nulls  are  assumed 
in  directions  Ur  75°,e=90° ) , (t>=105°,"=90°)  for  the  E°  field  and  U=60°, 
o=90°)  for  the  E1'  field. 

The  element  directions  for  the  first  array  are  <t>1=0»H1=(  i-1  )ti/36 
while  for  the  second  are  4>1=O,01  = n/6.  The  corresponding  optimum  polari- 
zation is  represented  in  Fig.  11. 


■ -v  -«•  • 


20 


o 

O 

0> 


Fig.  1.  Normalized  E°(<t>,n/2)  and  Ef(<j,  /2)  fields  of  a 1 inear 
array  with  6 elements  in  equal  distance  \/2  and 
direction  o’~0,  (Maximum  gain  subject 

to  the  constraint  on  nullsj. 
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V jV  i*  i.'f,' 


Fig.  2.  Normalized  E "($,v/2)  and  E^t,"/?)  fields  of  a linear 
array  with  10  elements  in  equal  distance  /2  and 
directions  0^=0,  /4.  (Maximum  gain  subject 

to  the  constraints  on  nulls,  sidelobe  levels 
and  efficiency  index.) 


2? 


array  with  6. elements  in  equal  distance  ^/2  and 


directions  o 1 = ( i - 1 ) / 1 0 . (Maximum  gain 

subject  to  the  constraint  on  one  null.) 
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| E*(<£,90*)|/|E®(90°f90°)|  AND  |E*l</>,90°)|/|E*(90o, 


4.  Normalized  E"($,n/2)  and  E*(4>,it/2)  fields  of  a linear  array 
6. element  in  equal  distance  \/2  and  directions  ij>1=0, 
o1  = (i-l  )tt/10.  (Maximum  gain  sub  jec  t to  the  constraints 
on  nulls,  sidelobe  levels  and  efficiency  index.) 


90  80  70  60  SO  40  30  20  10  0 


Normalized  E *(<t>,v/2)  and  E^(^,n /?)  fields  of  a linear 
array  with  6. elements  in  equal  distance  */2  and 
directions  $1=0,  o1=(i-l)  /10.  (Maximum  gain 
subject  to  the  constraint  on  nulls.) 
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Fi<).  6.  Normalized  E"($,n/2)  and  E (;,/2)  fields  of  a linear  array 
wjth  10  elements  in  equal  distance  74  and  directions  ^=0, 
’ = (i-l  )i/18.  (Maximum  gain  subject  to  the  constraints  on 
nulls,  sidelobe  levels  and  efficiency  index.) 
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|E®(<M0°)|/|E®(90°,90°)|  and  |E*(<£,90o)|/|E*(90o,90°)| 


r 


Fig.  7.  Normalized  E,J(<f>,7T/2)  and  E'*’((j>,n/2)  fields  of  a linear  array 
with  10  elements  in  equal  distance  x /4  and  directions  $i=0, 
01  (i -1  )tt/18.  (Maximum  gain  subject  to  the  constraints  on 
nulls,  sidelobe  levels  and  efficiency  index.) 


27 


1.8 


80° 


d/x 

(a) 


d/x 

(b) 


Fig.  9.  Optimum  polarization  versus  the  interelement 
distance  d/  tor  a 10-element  linear  array 

(a)  Constraint  nn  nulls 

(b)  Constraint  on  sidelobe  levels. 
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PHASE  PHASE 


d/x 


(a) 


(b) 

Fig.  10.  Maximum  directivity  versus  the  interelement 
distance  d/i  for  a )0-element. linear  array 
(Constraint  on  nulls),  (a)  <t>’=0, 
eMl-l)n/36.  (b)  e'=( i-1  )tr/6. 
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d/x 

(a) 


d/x 

(b) 

Fig.  11.  Optimum  polarization  versus  the  interelement  distance 
d/\  for  a 10-element  linear  array.  (Constraints  on 
nulls)  (a)  .j1  =0 , o1  = (i-l  )tt/36,  (b)  ^=0,  01  =(  i -1  )n/6. 
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VII. 


CONCLUSIONS 


Optimization  methods  have  been  presented  for  arrays  with  nonparallel 
wire  antennas.  It  has  been  shown  how  a performance  index  can  be  optimized 
subject  to  constraints  on  the  nulls,  on  the  sidelobe  levels  and  on  other 
indices.  It  has  also  been  shown  how  we  can  find  the  polarization  in  the 
direction  which  an  index  is  maximized.  Several  examples  were  included 
to  show  some  possible  uses. 


APPENDIX  I 


OPTIMIZATION  OF  ONE  PERFORMANCE  INDEX  SUBJECT 
TO  CONSTRAINT  ON  OTHER  INDICES 


Suppose  that  we  have  the  index  ,-j  which  must  be  maximized  subject 
to  constraints  on  other  indices  y2 »Y3 . • • * Index  y,  is  of  the  form 

[V]* C [ Bj  ] * [B^ ]+[Bg] * [B2]] [r] 

)]  = [^]*[M2][V]  84 

and  the  other  indices  are  given  as  follows: 

[!/']*[  Mi][V'] 

y,  - ~~  = a given  value  . (85) 

1 [V']*[Mi+1][V] 

The  conditions  Eg.  (85)  can  be  written  in  the  form 

[V'][MMi][V']  = 0 (86) 

where 

[MM-j  ] = Y-j  [Mi  +■]  ] - [Mi  ] . (87) 

For  maximizing  y^  we  will  use  the  Lagrange  method.  That  means  we  will 
find  the  maximum  of  the  function  L. 

[tf'3*[[813*[B1jHg2]*[B23][V] 

[V‘  ]*[M2][V  ] 

N 

♦ [ Ai[V'][MMi][V]  . (88) 

i =2  1 1 

Taking  the  first  variation  of  L,  XL,  egual  to  zero  we  have: 


<SL 


r . 1[[b1]*[bi3^23*[b2]]CV3 

r f s V I * "I  ^ _ 

[vmjtv] 


L*  J L'-'2  - 

[ vj  * [ [ Bj  ]*[  b 1Hii2  3 * [ B2  ] ] [ v ] 

~l7v']*[M2][v']|2 


[m2][v]  + 


N 

+ l A [MM.3CV3 

i=2  1 1 


+ {'’'}*  [6V‘]  = 0 


(89) 


where  { t*  signifies  the  conjugate  transpose  of  the  first  vector  { }. 
From  Eg.  (89)  we  have: 


[[fe1]*[B1]-^[B23*[B23][V] 

[V‘]*[M2][V7T~ 

[V‘3* [ [ B , ] ■ * [ B j3 +[ y * [ B 2 ] ] [ V1  ] 
I [V 1 ]*[m2][v ‘ ] 1 2 


[M2][V]  - 


N 

I *i [MM ■ ][V 1 ] 
i=2 


If  we  put 


u - ■= 


[b2][v] 

ra 


(90) 


1 

and  we  need  the  relative  values  of  the  voltages  of  Eg.  (90)  we  find: 


or 


[V1]  = 


[M2]  + l Pi [MM. 3 


-1 


[[B] ]*+u[B2]*] 


(91) 


[V]  = 


[m23  + pz[y2[m43-[m333  + ••••  -1  [[BilMB^]*]  . (92) 
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